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Recently the lower bounds of the intergalactic mag netic fields 10" 16 ~ I0" 20 Gauss are set by 
gamma-ray observations while it is unlikely to generate such large scale magnetic fields through 
astrophysical processes. It is known that large scale magnetic fields could be generated if there 
' exist cosmological vector mode perturbations in the primordial plasma. The vector mode, however, 

has only a decaying solution in General Relativity if the plasma consists of perfect fluids. In order 
, to investigate a possible mechanism of magnetogenesis in the primordial plasma, here we consider 

£\j ■ cosmological perturbations in the Einstein- Aether gravity model, in which the aether field can act 

as a new source of vector metric perturbations and thus of magnetic fields. We estimate the angu- 
lar power spectra of temperature and B-mode polarization of the Cosmic Microwave Background 
(CMB) Anisotropies in this model and put a rough constraint on the aether field parameters from 
latest observations. We then estimate the power spectrum of associated magnetic fields around the 
recombination epoch within this limit. It is found that the spectrum has a characteristic peak at 
k — 0.1/iMpc -1 , and at that scale the amplitude can be as large as B ~ 10~ 22 Gauss where the 
upper bound comes from CMB temperature anisotropies. The magnetic fields with this amplitude 
can be seeds of large scale magnetic fields observed today if the sufficient dynamo mechanism takes 
place. Analytic interpretation for the power spectra is also given. 
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PACS numbers: 04.50.Kd, 98.80.Es 



I. INTRODUCTION 



Astronomical observations have shown that magnetic fields exist ubiquitous in various astrophysical objects, ranging 
— i ' from planets to clusters of galaxies or even larger systems [U-Il] . On large scales such as in galaxies and clusters of 
t> . galaxies, measurements of the Faraday rotation indicate that the strength of magnetic fields is about ~ 1CP 6 Gauss. 
Recently the lower bounds of the intergalactic magnetic field are found as 10 -16 ~ 10~ 20 Gauss [1, Hj], depending on 
the data used and assumptions they have made. 

The origin of such large scale magnetic fields remains as an enigma and much effort has been made to explain the 
origin in a wide variety of ways. It is believed that the magnetic fields in galaxies can be amplified and maintained by 
the dynamo mechanism, which is a hydro-magnetic process with magnetic reconnections (e.g. refs. [1, However, 
we still need a seed field for the dynamo process to act on. To explain the observed magnetic fields ingalaxies in 
the present universe, the seed field should be as large as 10~ 20 ~ 10~ 30 Gauss at kpc comoving scale Q- It is an 
interesting scenario that such small seed fields may end up with the intergalactic magnetic fields of > 10~ 20 Gauss 
due to adiabatic compression. 

The seed magnetic fields are possible to originate from quantum fluctuations of the electro-magnetic fields stretched 
by inflation in the early universe. If the conformal invariance is broken by some mechanisms during the inflation era 
magnetic fields with a large coherence length can be naturally generated beyond the horizon scale. In this 
case, however, there are associated nagging problems, namely, the back reaction and the strong coupling problems. 
The former is that inflation fails to proceed if the generated electromagnetic fields dominate the energy density in 
the universe during inflation. The latter problem is related with the naturalness of the model building. Due to 
these problems, there is no satisfactory model so far to explain observed magnetic field s fl2l. Nevertheless, the 
effects of primordial magnetic fields on the present observations have been well investigated [12h15| | . Another possible 
mechanism to generate seed magnetic fields is the phase transition in the early universe (e.g. ref. [6[). Phase transition 
releases the free energy as a latent heat with forming bubbles, and the bubble collisions can generate electric current 
and hence magnetic fields. However, the magnetic fields generated at the phase transition generally have small 
coherence length which corresponds to the Hubble radius at the transition. In this case, therefore, inverse cascading 
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processes of magnetic fields are necessary to explain the ones at large scales. 

On the other hand, some astrophysical processes may become natural candidates of the origin of the seed magnetic 
fields. For example, the Biermann battery [1(1], which is effective in non-adiabatic situations such as shocks, can 
generate magnetic fields in starts [I?} , supernova shocks [3, [l9| , large scale structure formation [2(| , and cosmological 
re-ionizatoin [2lj ]. However, since the mechanism works effectively only in places where the matter or astronomical 
objects are present, it may be difficult to explain the existence of the intergalactic magnetic fields if they are in void 
regions (22j. 

In this paper, we consider yet another mechanism to generate the seed magnetic fields, i.e., the vorticity in the 
primordial plasma before the recombination epoch [23| . Because photons push electrons more frequently than ions 
through Thomson scatterings, the vorticity in the photon fluid can induce circular current and thus generate magnetic 
fields. The problem here is that in a Friedmann universe with perfect fluids, the solution of the vorticity has only a 
decaying mode. In order to have non-zero vorticity, some mechanisms are proposed, which include the free-streaming 
neutrinos 0, [H[ , and the non-linear couplings of the first order density perturbations (25l - [27| . 

Following this idea, here we show that a possible modification of gravity, namely the Einstein- Aether gravity [28| , can 
also generate non-zero vorticity and thus magnetic fields. The Einstein- Aether gravity is known as a healthy extension 
of Hof ava-Lifshits gravity at low energies [2^, H3| , which originally proposed by Hof ava [3l|, [32j as a candidate of the 
theory of quantum gravity. The Einstein- Aether gravity contains a new regular vector degree of freedom which is 
called "the aether field" . Effects of the aether field have been discussed intensively, for example, in connection with the 
current accelera ting exp ansion of the univers e I33h35| . Cosmic Microwave Background radiation (CMB) temperature 
anisotropies [HI l36l l37l| and compact objects [38 -40] . The purpose of this paper is to examine the role of the vorticity 
excited by the aether field in the generation of magnetic fields before the recombination epoch. 

The plan of this paper is as follows. In section II, we review the Einstein- Aether gravity and describe a formalism of 
the vector-mode perturbations. In section III, we focus on the evolution equation of magnetic fields in the primordial 
plasma. In section IV, we explore the evolution of the aether field and perturbations. Our main results will be 
described in this section. Section V is devoted to our summary. In Appendix A, we summarize the observational 
and theoretical constraints on the aether parameters. In Appendix B, we formulate cosmological perturbations in the 
Einstein- Aether gravity. 

Throughout this paper, we use the units in which c = K = 1 and the metric signature as (— , +, +, +). We obey 
the rule that the subscripts and superscripts of the Greek characters and alphabets run from to 3 and from 1 to 3, 
respectively. 

II. EINSTEIN-AETHER GRAVITY 

In this section, we summarize the Einstein-Aether gravity. The action for the Einstein-Aether gravity is given 
byHI 

S = j d^xj—g [R + C A }+ J d 4 x^C m , (I) 

C A =K + \(A a A a + l) , (2) 
K = K<*% (V a A») {VpA") , (3) 

K a %, = drt^ + c 2 5^Si + c 3 S a J^ - c A A a A' 3 g^ . (4) 

Here, C m is the Lagrangian of the ordinary matter, R is the Ricci scalar, and Ca is the Lagrangian of the aether 
field A^. We assume that the aether field does not couple to the matter. The constant G is different from Newton's 
gravitational constant which is locally measured. That is to say, G is a "bare" parameter. We will see that G is 
renormalized. The coefficients C\, C2, C3 and C4 are the set of parameters of the Einstein- Aether gravity. These 
parameters are constrained by observations and theoretical hypotheses as summarized in Appendix A. This theory 
contains a massive ghost [ill ], and one usually imposes a fixed-norm constraint on the vector to remove such a ghost. 
This is equivalent to forcing a Lorentz-breaking on the vacuum expectation value. The parameter A is a Lagrange- 
multiplier to fix the norm. 



where Ca is written in the form: 
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To derive a set of equations of motions, we take variations with respect to variables. Variation with respect to the 
Lagrange-multiplier A imposes the condition that 

A"A„ = -1 . (5) 

This equation is the constraint equation of the aether field . Owing to the fixed norm, we will eliminate the 
possibility that the degree of residual freedom acts as a ghost [42| . 

Variation with respect to the aether field gives the equation of motion and the continuous equation for the 
aether field as 



V Q J a + Ci A a (V Q A^) (V M A 7 ) = XA^ , (6) 



where we define J a ,, as 



J\ = K a %V p A» 

= Cl V a A^ + c^WpA^ + c 3 V^A a - c 4 A a A^V^ . 

Variation with respect to the metric g^ v gives the Einstein equation with the aether field, 



(7) 



R„u - \R9»v = T$ + SttGTW , (8) 



where we also define energy momentum tensors Tffl and T^f^ as 

rp(A) = 1 S(^jC A ) 



T (M) 



2 8b/=gC 



M I 



(9) 



Note that the factors of the energy momentum tensors for the aether field and the matter arc different. The concrete 
expression of T^t can be written as 

^ = V. (2%) + y& - 4t 4) + ^A v + \g^L A . (10) 

Here we define the tensors X" , yffi and as 



Zj& = c i A a A 13 (V Q A M ) (y p A v ) , 

yjg) = ci (yApVvAv - v^vm , (n) 



"iiv = A^J^ — A^J 17 u j — A a J^v) 



where the parentheses, , mean symmetrization. For convenience we use the abbreviations, which are written in 
the forms: 

a = a + 3c 2 + c 3 , 
C123 = ci + c 2 + c 3 , 
C13 = Ci + C3 , 
Cl4 = Cl + C 4 . 

Because theoretical parameters have four independent components, we will treat a, ci 4 , C13 and ci as the independent 
model parameters. 



A. Background cosmology 

In this subsection, we explore background cosmology where the space-time is given by a spatially flat Friedmann- 
Robertson- Walker metric. The line element is given by 

ds 2 =a 2 (r/)[-(i?7 2 +dx 2 } . (13) 
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From Eq. ([5]), the background components of the vector are given by 

= (a" 1 , 0,0,0) . (14) 
Substituting Eq. (fT4")l into Eq. ^ , the p = component of Eq. © can be written as 

3 



A = 

a 



s (ci 23 n 2 - c 2 n) , (15) 



where H = a/ a, and a dot represents a derivative with respect to the conformal time rj. Similarly we substitute Eqs. 
(IT31 and (fl4"|) into Eq. ([5]), and we derive the Friedmann equations with the aether field as 

n 2 = ^^a 2 p , 

Here we define G cos = x G a j 2 i an d p and p are the energy density and pressure of the ordinary matter, respectively. The 
effect of the aether field can only be seen in the renormalization of Newton's gravitational constant in the background 
cosmology. 

In the Newtonian limit which should be applied to laboratory measurements, the Einstein equation can be rewritten 
in the form of the Poisson equation [43j as 

V 2 $ = 47tGnP™ , (17) 

where Gn = , t G to - Because G cos and Gn can be different, the ratio of G cos and Gn is constrained by observations 
such as the light element abundance from the Big-Bang Nucleosynthesis (BBN) [43j (see Appendix A). From that we 
have a constraint on the aether parameters as C14 + a < 0.2. We assume that C14 = —a in this paper for simplicity, 
which leads to G cos = Gn- 

B. Vector- mode perturbations 

In this subsection, we consider cosmological perturbations in the Einstein- Aether gravity (see also Appendix B). 
Since we are interested in generation of magnetic fields via the vector-mode perturbations, we focus on the vector- 
mode in this paper, and we omit indices ^ ±1 - ) or ( v ' that represent the vector-mode in the Appendix. B. We will work 
in the synchronous gauge in which the metric is given by 

ds 2 = a 2 [-drj 2 + (Sij + hij) dx i dx j ] , (18) 

where is the metric perturbation. Under this metric, the aether field can be expressed in the form as 

'1 V 



A^ 

' a a J ' (19) 

A^ = {-a,aVi) , 

where we have imposed A^A^ = — 1. Raising or lowering indices of Vi is done by 8^. Substituting the above equations 
into Eqs. ([8]) and ©, we derive the evolution equations of the perturbed quantities. The equation of motion for the 
aether field is given by 

C14 (V + 2HV) - {a - c 14 )HV +(a + c 14 )H 2 V + ak 2 V + ^c 13 k 2 a = . (20) 
The Einstein equations are given by 

& + 2Ha = - j-^ — (v + 2HV) - , (21) 

1 + C13 V / (1 + Ci3jk 

k 2 G = \l6irGa 2 q - c 13 k 2 V] , (22) 

1 + C13 
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where a = h/k denotes the shear, and pn and q denote the anisotropic stress and the heat flux of the ordinary matter, 
respectively. 

To see the behavior of the aether field, we solve Eq. (|20[) without the ordinary matter contribution on a first. The 
ordinary matter is subdominant in the evolution of the vector perturbations in the Einstein- Aether gravity, although 
we will include the matter contributions correctly in our numerical calculations. This is because the ordinary matter 
does not couple to the aether field while the aether field and the metric couple directly. Thus the aether field can 
amplify predominantly the metric perturbation as its source, namely, Eq. (|22[) can be rewritten as 



ct~ ^—V. (23) 

I + C13 

This assumption is justified by the initial conditions and the numerical results. Using Eqs. (|20l) and (I23|) . we can 
obtain the evolution equation of the aether field in the closed form as 



V + 2HV 



1 + ^ )U 2 + (l-±)H 

CuJ V C14 



V + 4k A V = , (24) 



where c v is the (effective) sound speed of aether field perturbation and defined by 

If 1 



c?, = — U-±_J2_j . (25) 

d4 V 2 1 + c 13 



As a consequence of the existence of the sound speed, the aether field will have a characteristic scale i.e., the "sound 
horizon" , which is defined as fcgH = aH/c v . Depending on whether the mode is inside or outside the sound horizon, 
the behavior of the aether field will change. 

When we assume a long-wavelength limit c v krj <C 1 in the radiation dominated era (RD) and the matter dominated 
era (MD), the solutions are given by 

F M °J!fT d R ° (26) 
1 {krj) v ™* MD , y ' 



where ^ ra d and ^ ma t are defined by 



_ -1 + v/1 - 8a/ci4 

t'rad = 7, 



_ -3 + v/1 - 24a/ci4 

"mat = ~ 



(27) 



These results show that the aether field has the power-law dependence on kr\ outside the sound horizon. Because we 
are interested in a non-decaying regular mode in the radiation dominated era, we assume that ajc\^ < 0. Taking 
into account the parameter constraint from Eq. (|A8|) . the power-law indices, z/ ra( j and v ma t, must be satisfied with the 
conditions < i-'rad < 1 and — 1 < ^ ma t < 1. 

Next, we assume the short-wavelength limit c v kr\ 1. In this limit, the solutions are given by 

vik^A^'^l R ° (28) 

V ' " |e lc " fc "/(fcr/) 2 MD . V ' 

Both in the radiation and matter dominated eras, the aether field decays as a a -1 with oscillations. 



III. MAGNETIC FIELDS GENERATION 

In this section, let us consider the generation of magnetic fields. We follow Ref. [25[ to formulate the generation 
process. When there exists the difference of velocities between baryons and photons, magnetic fields can be generated 
[25l - l27l . |44l |45j |. This is caused by the difference of the Thomson cross sections between electrons and protons. In 
other words photons push electrons more frequently than protons. As a consequence, the charge separation appears. 
Photon's bulk pressure creates the charge separation which induces electric fields, and then magnetic fields will be 
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generated as well. The evolution equation of magnetic fields is derived by the combination of the Maxwell equations 
and Euler equations for electrons and protons with Thomson scattering [25j. It is written in the form: 

im.lfI|2£ e .^ K) _ %) , ( 29) 

where <tt is the Thomson cross section and e l ^ k is the Levi-Civita symbol and p 1 is the energy density of photons 
and i> 7 or Vb are the velocity of photons or baryons. When we consider scalar perturbations for the velocity fields, 
the right-hand side of the above equation should vanish. In this paper, we solve the above equation with an initial 
condition B = at z = 10 9 , which roughly corresponds to the time of neutrino decoupling. By integrating the above 
equation, the square of the magnetic fields is given by 

a 4 £ l (M)B*(k',t) = (^pj (S^S 1 ™ -6 jm 5 M )k k k' m J a(t')p 1 {t')6v J (k,t')a{t")p. ( (t")5v* e (k',t")dt'dt" , (30) 

where 8vj(k,t') = w 7 j(k, t') — Vbjik, t'). By taking an ensemble average and defining the power spectrum, we obtain 

2?r 2 

(^(k,i')^(k',t")) = (27r) 3 —Vv(k)P J e(k)5v(k,t')5v(k,t")S 3 (k~k r ) , (31) 

P je (k) = S jt - kjk e , (32) 
where Sv = u 7 — Vb is the transfer function of the baryon-photon slip, and Vv{k) is defined by 

(v< A >fa,k)V*Wfa,k')) = (2n) 3 ^V ( v X \k)6(k-k') , (33) 

where A means two helicity states. Because the parity is not violating in the Einstein-Aether gravity, two helicity 
states A = ±1 should not mix. Thus we omit the superscript ^ hereafter. The power spectrum of magnetic fields is 
defined as 

( J B i (k,t)S*(k',t)) = ^SBik^t^ik-k') . (34) 
Then we can express SB{k,t) by the initial power spectrum for the aether field T-V(fc) as 



A- 3 



l (t)^S B (k,t)=f^f]\vv(k)k 2 



dt a(t )p-y(t )Sv(k, t') 



2 



(35) 



The initial power spectrum may depend on the inflation model considered. Generally, the initial power spectrum is 
assumed to be given by a power law as 

V v {k)=A v " , (36) 

where k — 0.002Mpc _1 . If we fix the inflation model and the subsequent reheating model, the initial amplitude Ay 
and spectral index n v can be expressed as a function of the aether parameters and inflation model parameters (see 
Appendix B). 

The amplitude of velocity difference between baryons and photons to generate magnetic fields is dependent on the 
amplitude of the metric perturbation of the vector-mode. The aether field induces the metric perturbation of the 
vector-mode dominantly, and therefore the evolution of the aether field is closely related to that of the magnetic fields 
and their power spectrum. 

To understand the behavior of the vector-mode velocities in the presence of aether field, we should solve the Euler 
equation for baryons, which is given by 

ib + U-Vb = --j^-anecrrivb - v~) , (37) 
3pb 

where n e is the electron number density. We perform a multipole expansion of the Boltzmann equation for photons 

as 

w 7 + —kir-y — —an e aT(v-y — Vb) , 
8 

8 , r S , ( 9 9 r, \ 8 , 

tt 7 + -kl 3 - -fcw 7 = ~an e a T I j^tTj - -E 2 I - -ka , (38) 
I f£ + 2 \ 

+ 21 + 1 1 ITi Ie+1 ~ ) = ~ aUe<7Tli ' 
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where Ii is the £th order moment of the photon distribution function and Eg is the £th order moment of E-modc 
polarization. The aether field does not affect the Boltzmann equation for photons and the Euler equation for baryons. 
However the aether field may change the initial conditions for the matter part because the Einstein equation contains 
the aether field. 

Before moving to our numerical calculations, we summarize the initial conditions and the set of equations under 
the tight coupling approximation. The initial conditions with the aether field have been derived by assuming that the 
universe is deep in the radiation dominated era and expanding the equations in powers of kr\ up to the lowest order 
46]. Because we arc interested in the effects of the aether field, we ignore the regular vector-mode in the presence of 
the neutrino anisotropic stress which is investigated in (25|. Then the aether field in powers of krj up to the lowest 
order is given by 

V = A kV ^ [l - (l - «f) f ] , (39) 

where uj — QmHo I V^Tr and &>i are ordinary cosmological density parameters. The initial conditions for the other 
variables are given by 

v-f = , 



w„ =0 



a = — - 



^rad C 13 
"r*ad + AR t 1 + C 13 



AkTf** , (40) 



Hi = 8 ^rad C13 , , 1+^ 

Pu 15(1 + v xaA ) i/ r * ad + m* u i + Cl3 

where R* = (1 — a/2)/(l + Cis)R v , R v = Q u /Cl r , u* ad = 5(1 + ^ ra d)(2 + f ra d)/2. We apply these initial conditions to 
our numerical calculations. 

Deep in the radiation dominated era, photons and baryons frequently interact with each other. These fluids are 
tightly coupled because the opacity f = an e UT is large. Hence the tight-coupling parameter e allows us to expand 
the equations. The tight coupling parameter can be expressed as 

where is the baryon density normalized by the critical density and h = i/o/100 (km s _1 Mpc -1 ) is the normalized 
Hubble parameter with Hq being the Hubble constant. We expand the equations using the tight coupling parameter 
up to second order. Then we derive the velocity difference between baryons and photons as 

k RH 4 fk\ 2 R 



Sv = — v y (v^ + cr), (42) 

r(l + R)k 7 15 \t J 1 + iT 7 ; ' v ; 

where R = |^ and we have neglected the cosmological redshift terms. In Eq. P2]l . although the first term does 

not depend on the wavenumber, the second term depends on k 2 . Because the initial conditions of our numerical 
calculations imply that a 3> w 7 , we can see that the second term dominates over the first term. In contrast, if we 
consider the vector-mode without the aether field, i> 7 and a are in the same order. In this case, the first term in 
Eq. (|42[) must dominate in the early universe [25]. This fact appears in the difference of k dependence of 5v as shown 
later. Namely, in the case without the aether field, we have asymptotic scaling as Sv cx fc°. In contrast, in the case 
with the aether field, the scaling relation becomes Sv cx fc 2 . 

IV. NUMERICAL RESULTS 
A. Time evolution 



In this section, we show the result of our numerical calculations. To solve the set of equations, namely, the evolution 
equation of the aether field (Eq. ([6])), the Einstein equation with the aether field (Eq. ©) and the Boltzmann equation 
for the photon distribution function (Eq. (|3"5|) ). we modified the CAMB code [471 ] . 
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In order to see the evolution of perturbed quantities, we show two results of each parameter set, and the magnitude 
of these parameter sets is extremely different from each other. First, we fix the aether parameters as 

C14 = -1.0 x icr 4 , 

c 13 - -2.0 x 10~ 4 , 
a = -c u = 1-0 x 10~ 4 , (43) 
Ci = -3.0 x 10~ 4 , 
c 2 v = 3.0 . 

These values are consistent with all the constraints coming from the observations and the theoretical hypotheses 
as shown in Appendix A. Time evolutions of perturbed variables are depicted in Fig. [TJ As mentioned earlier, the 




FIG. 1: Time evolutions of perturbed variables in the Einstein-Aether gravity for the aether parameters fixed as Eq. (|43|) at 
k = 10 /tMpc -1 . Shown in the left panel are the shear a, the aether field V, the velocity of baryons Vb, the velocity of photons 
v 1 and the velocity difference between baryons and photons Sv. Black solid vertical line does not express a horizon entry time 
but "a sound horizon entry time" . The metric perturbation of the vector- mode a and the aether field start to decay not after 
the horizon crossing but after the sound horizon crossing. Subsequently, they start to oscillate and damp. In the right panel 
time evolutions of the velocity difference between baryons and photons Sv are shown for some wavemimbers as indicated. We 
normalized by the factor of k 2 to remove the scale dependence in early time (see the last paragraph of the previous section). 



perturbation of the aether field V and the metric perturbation a dominate over the other perturbation variables in the 
early universe. When the modes of Sv and a enter the sound horizon, they start to oscillate and decay in proportion 
to the inverse of the scale factor. On the contrary the baryon and the photon velocities increase until breaking the 
tight coupling. 

When the Fourier modes reach the Silk damping scale, perturbations in the photon fluid start to vanish exponentially 
along with the baryon velocity and the potential a in the absence of the aether field (25j . In contrast, in the presence 
of the aether field, such exponential damping does not appear as shown in Fig. [1] This is because the aether field 
can keep the vector potential a large enough, and push the photon fluid in equilibrium with a through the photon 
anisotropic stress even in the photon diffusion regime. The baryon velocity decays faster than the photon velocity in 
a particular case shown in Fig. [T] (green dashed line). This can be understood by noting the fact that baryons are 
non-relativistic having no anisotropic stress, and therefore their velocity is completely determined by the Compton 
drag force from the photon fluid independently of the metric perturbation a (see Eq. (|37p). Because the photon 
velocity is sustained by the aether field and kept large enough even after the tight coupling is broken down, the 
equation for the baryon velocity Eq. (|3"T)) is reduced to 

Vh « —^-an e GTV 1 . (44) 

Because the right hand side of the above equation starts to damp exponentially with rapid oscillations right before 
the recombination epoch, the baryon velocity also damps exponentially with oscillations. As a result, the difference 
of the velocities between baryons and photons has the same amplitude as the photon velocity at late time. 

In Fig. [I] before crossing the sound horizon, the evolution of Sv may be explained by use of the tight coupling 
approximations, i.e., Eqs. ([22]), (|41~j) and (|42|) . Since a and the aether field dominate the perturbation variables, we 
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can ignore the photon velocity in Eq. (I42[) . Then we have an approximate expression as 



8v 



A_ fk 
15 \f 



R 



1 + R 



(45) 



Substituting the time dependence of f , R and a during the radiation dominated era, the above equation is reduced 
to the simple form as Sv oc k 2 a 5 a oc fc 2 a 5+lyrad . After crossing the sound horizon, Sv continues to grow until the tight 
coupling broken. Subsequently, the velocity difference between baryons and photons Sv decays together with a as 
discussed above. 

To see the dependence of the aether parameters, we consider another parameter set: 



ci4 = -0.9 x 1(T 1U 
ci 3 = -1.8 x 1(T 10 
a = -cm = 0.9 x 10~ 10 , 

ci = -0.9 x 10~ 10 

cl = 1.0 . 



(46) 



The result is depicted in Fig. [T] and Fig. [2] The relation between V and a obtained from Eq. 



is indeed satisfied. 
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FIG. 2: Same as Fig. [T] but for the different aether parameters fixed as Eq. (|46 



We find that 6v in Fig. [5] is smaller than one in Fig. [T] If C13 <C 1, we can rewrite the equation ([23]) as a ~ —C13V. 
Because of this relation, a is suppressed by the factor of C13 compared with V. Furthermore, we can see Sv ~ u 7 ~ cr 
at late time. Therefore, magnetic fields are also suppressed at late time by the factor of C13 if C13 <C 1. 



B. Power spectra of temperature and B-mode polarization 



In this subsection, we calculate the power spectra of CMB temperature and B-mode polarization. We assume that 
the initial power spectrum of the aether field is given by a power law as Eq. (|36p. Although the initial amplitude Ay 
and spectral index n v may depend on the early universe model, following the discussion in Ref. [46l |. we assume that 

*v = 3 - Jl - — 4 % 2 , (47) 

where e = l — H'/H 2 is the slow-roll parameter. In addition, we assume that the spectral index does not change in the 
scales of interest. The combination of the aether parameters a/cn is constrained as — 1 < a/cn < for perturbation 
of the aether field to have a growing mode solution in the radiation dominated era and the isocurvature mode does 
not grow, see Eq. (|A8|) . The spectral index is also constrained as 2 — 2e/(l — e) 2 < n v < 2 for e <C 1, and thus we 
fix the spectral index as n v = 2 in our numerical calculations for simplicity. Since the aether amplitude may further 
change during the reheating stage, we treat the initial amplitude as a free parameter (irij . 
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Before moving to the calculation of the magnetic field spectrum, we can give a rough constraint on the initial 
amplitude and the aether parameters from the CMB temperature anisotropics. The CMB temperature (TT) and 
B-mode polarization (BB) power spectra for the fiducial parameter set are depicted in Fig. [3J From the figure we find 
that the TT power spectrum from the aether field looks similar to that from the primary tensor perturbations, with a 
slightly different horizon scale, i.e., the first peak location, as discussed in section IV. A. Therefore the constraint on 
the aether field from the TT power spectrum comes from the low multipole components, which is similar to the case 
of the primary tensor modes. As for the parameter dependence of the TT and BB power spectra, we find that the 
amplitudes of the spectra depend on the combination of the initial amplitude Ay, the aether parameters C13 and c v 
as Ayc\ 3 c~ 4 . One can understand this dependence as follows. First of all, the amplitude of the metric perturbation 
a is proportional to C13 as shown in Eq. (j2"3")l . and thus Cg oc a 2 oc c 2 3 . Next, from Eqs. and ([25]). the aether 

perturbation V evolves as V oc r\ outside the sound horizon during the radiation dominated era while V decays 
as V oc ?7 _1 with oscillations inside the sound horizon. The resultant amplitude inside the sound horizon can be 
expressed as V — v\ n i (?7*/77ini) (v/v*) 1 where the subscript "ini" means the initial value and 77* represents the epoch 
of sound horizon crossing, which is given by 77* = c" 1 . Therefore V oc rjl = c^ 2 and the amplitude of the spectra 



will be changed in proportion to c v 4 . 
dependence explicitly in Fig. @] 



Note that the peak positions of the spectra also depend on c v . We show this 
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FIG. 3: The CMB temperature and B-mode polarization anisotropy power spectra. "Primary Scalar" and "Primary Tensor" 
indicate the primary spectra from the scalar and tensor perturbations in the standard cosmology. "Aether Vector" indicates the 
spectra from the vector mode perturbation in the Einstein-Aether gravity. The aether parameters are set as C14 = —1.0 x 10~ 4 , 
C13 = —2.0 x 10 -4 , a = 1.0 x 10 -4 , ci = —3.0 x 10~ 2 and c 2 , = 3.0. The initial amplitude and the spectral index are fixed as 
■4yCi3C^" 4 = 2.0 x 10 -17 , and n v = 2.0 respectively. We assume that the tensor-to-scalar ratio is r — 0.1. 



From the above discussion and Fig. [3J we find a new constraint on the aether parameters and the initial amplitude. 
Because the dependence of Cf ethcl ' on the aether parameters is proportional to Ayc\ 3 c~ A , we find a relation from 
Fig.[3jas 



C aether 
t TT : 



' tt 



Ayc 2 13 c v 4 s 

1-17 



v 2.0 x 10- 



(48) 



where Cf^lf is the TT power spectrum by primordial gravitational waves with the scalar-tensor ratio r = 0.1. 



By imposing C : 



act lie 



<c 



GW(r=0.1) 
TT 



which is the current upper bound obtained by the WMAP team [48|, |49( , we can 



place a constraint on the aether parameters and the initial amplitude as 



-4yc 2 3 c- 4 < 2.0 x 10~ 17 . 



(49) 



When we assume the single field slow-roll inflation, the initial amplitude Ay is given in the aether parameters (See 
Appendix B). If we adopt Eq. (|B32j) . the above inequality can be rewritten as c 2 3 (— ci4)~ 1 c~ 5 < 10 43 . 

It should be noted that the B-mode power spectrum shows distinctive feature of oscillations at i > 80, and thus 
future precise measurements of B-mode polarizations can be used to discriminate the contributions from the aether 
field and the primary tensor perturbations (Ref. [46] ) . Note also that the aether field, with this upper bound on 
Avcl 3 c~ 4 , gives negligible contributions to the TE (E-mode and temperature cross correlation) and EE (E-mode 
auto correlation) power spectra compared to those from the standard (observed) density perturbations, and we have 
omitted them here. 
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FIG. 4: The CMB temperature anisotropy (TT) power spectra with different aether parameters; C13 and c„. Here we impose 
the condition, Avc^c^ 4 = 2.0 x 1CP 17 , on the aether parameters and the initial amplitude. "Primary Tensor" indicate the 
primary spectrum from the tensor perturbation in the standard cosmology. The spectra from the aether perturbation are 
depicted by solid lines, (ci 3 , c 2 v ) = (-2.0 x 10" 4 , 3.0), (-2.0 x 10" 6 , 9.0) and (-4.0 x 10" 10 , 27.0), respectively. This figure 
suggests that if AvCi3C^ 4 is fixed, the amplitudes of TT power spectra are almost the same at low multipoles. On the contrary, 
the peak locations are shifted. The TT power spectrum from tensor modes with r — 0.1 is also plotted as a reference (broken 
line). 




k[hMpc l ] 



FIG. 5: The spectrum of magnetic fields with c 14 = -1.0 x 10~ 4 , C13 = -2.0 x 10~ 4 , a = 1.0 x 10~ 4 , ci = -3.0 x 10 -4 
and <? v = 3.0 at 1 + z — 1100, 1 + z = 2000 and 1 + z — 3990. The initial amplitude and the aether parameters are fixed as 
Avci 3 c~ 4 = 2.0 x 10" 17 . 



C. Power spectrum of magnetic fields 

The final results of generated magnetic fields spectra are depicted in Fig. [5] This figure shows that there are three 
characteristic scales. Let us focus on the scales, 10~ 4 < k/hMpc^ 1 < 10 -2 , which are outside the sound horizon at 
1 + z = 2000, right before the recombination epoch. In these scales, the magnetic field spectrum is proportional to 
k A . Next, let us focus on the scales which are inside the sound horizon and in the tight coupling regime (Sv <C w 7 ) 
, 10~ 2 < fc/ZiMpc" 1 < 10°. In these scales the aether field begins to oscillate after the sound horizon entry. The 
velocity difference Sv, the source of magnetic fields, also begins to oscillate because Sv is dragged by the aether field 
through a (see Eq. (|45|) . The velocity difference grows with oscillations until the tight coupling breaks down. In these 
scale, the spectrum of magnetic fields is proportional to k 1 with oscillations. 

Finally we focus on the scales, k/hMpc^ 1 > 10°, in which the tight coupling approximation is no longer valid. In 
these scales Sv ~ i> 7 ~ a, and Sv begins to decay as cx a -1 together with a. Consequently, the spectrum of magnetic 
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fields decays with oscillations at k/hMpc^ 1 > 10°. As a result the spectrum of magnetic fields becomes proportional 
to fc _1 with oscillations. 

At l+z — 1100, the amplitude of magnetic fields has an extra enhancement due to the recombination process. When 
the Fourier mode is inside the sound horizon and in the tight-coupling regime, baryons oscillate with photons. The 
Fourier modes in this scale can have an extra growth because the tight coupling becomes weaker as the recombination 
process takes place. On the other hand, there is little enhancement for the Fourier modes at scales where the tight 
coupling has already been broken down at recombination. As a result, around the recombination epoch, the position 
of the peak of the spectrum moves slightly toward the smaller wavenumber. 




FIG. 6: The spectrum of magnetic fields with difference of the aether parameters and the initial amplitude at 1 + z = 1100. 
To normalize the amplitude of the spectra, we multiply by <? v . We impose only the condition, this condition is same as Fig. [4] 
AvCi^c~ 4 = 2.0 x 10~ 17 , on the aether parameters and the initial amplitude. The spectra of the magnetic fields are depicted 
by solid lines, (ci 3 , c\) = (-2.0 x 10~ 4 , 3.0), (-2.0 x 10~ 6 , 9.0) and (-4.0 x 10 -10 , 27.0), respectively. 

Let us now estimate the largest amplitude of the magnetic fields allowed from the current observations. We find that 
the amplitude of the generated magnetic fields also depends on the combination of Avcl 3 c~ a . This aether parameter 
dependence can be understood in the same way as CMB power spectra as follows. Here we assume the universe is 
dominated by radiations. Magnetic fields are given as the time integral of the velocity difference Sv (Eq. (pJ5"|) ) as 

^ a*(t)-^S B (k,t) = (^) V^Mkjk J" dt'a(t')p 7 (t')Sv(k,t') 

cx [ V dr)'r)' 2 p^(r)')5v(k,r)') (50) 



JO 

cx / drj 'if Sv(k, rf) . 
Jo 

Time dependence of Sv has two characteristic epochs; the sound horizon crossing and the tight coupling breakdown. 
Before the sound horizon crossing, Sv evolves as Sv cx rf as shown in Figs. [T] and [2J (also see Eq. (05])). Once 
the Fourier mode enters the sound horizon, Sv evolves as Sv cx rf with oscillations. When the tight coupling breaks 
down, moreover, Sv evolves as Sv cx r\~ x with oscillations. Therefore, by substituting the above time dependence 
into Eq. (|50)) . we obtain the time dependence of magnetic fields as B cx rj 5 , rf jc v , and ij~ 3 /c v , before and after 
sound horizon crossing, and after tight coupling breakdown, respectively. Note that the above dependence of rf /c v 
and rj~ 3 /c v are obtained from the fact that / drjrf 1 e™ CvT] ~ rf 1 j (ikc v )e lkc " J ' n ~ rf L /c v . Accordingly, the dependence 
of magnetic fields on the aether parameters can be understood in the same way as CMB power spectra. B cx 
Bi n i(r]* / r] ini ) 5 ((?7tc/??*) 2 /c v ) ((r] / Vtc)^ 3 / c v ) cx c~ h c~ x (? v cx c~ 4 where 77* cx c" 1 and ry tc oc c" 1 are the epochs of sound 
horizon crossing and tight coupling breakdown, respectively. Since the amplitude of Sv is proportional to C13, the 
overall dependence is given by B 2 cx Avc( 3 c~ s , with Ay being the initial power spectrum amplitude. This dependence 
is explicitly depicted in Fig. [5] If Ayc\ 3 c~ % is fixed, the amplitudes of generated magnetic fields are almost the same. 
Note that, however, peak positions of the magnetic field spectrum depend on c v as the CMB temperature anisotropy 
power spectrum. 

Let us estimate the maximum amount of magnetic fields allowed from the limits of the CMB power spectra Eq. P^)l 
and the constraint on the sound speed c v (Appendix A. 5). We find that the spectrum of magnetic fields with aether 
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field parameter dependence is given by 

4 

[Gauss] (fc/frMpc- 1 < lCT 2 /^,) 



[Gauss] (lO~ 2 /c v < k/hMpc' 1 < l/c v ) (51) 

[Gauss] (1/c < fc/ZiMpc- 1 ) . 

If we set AvCi 3 c~ 4 = 2.0 x 10 -17 (Eq. (fifTJl) ]) and c v — 1.0, which satisfies the constraints in Appendix A, magnetic 
fields have the largest amplitude as 10~ 22 Gauss at 0.1 ftMpc -1 . 
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V. SUMMARY 



In this paper, we explored a mechanism of the generation of magnetic fields in the Einstein- Aether gravity. This 
theory contains a dynamical vector field, i.e., the aether field, and the aether field excites the metric perturbation of 
cosmological vector modes. Because the metric perturbation generates the velocity difference between baryons and 
photons, magnetic fields are naturally generated through Thomson scattering. 

We investigated evolution of vector perturbations in detail in the Einstein Aether gravity. We found that the 
evolution of the vector metric perturbation a is quite different from the one with vector perturbations in General 
Relativity 25]. In particular, in the case with the aether field the Silk damping does not arise even after the tight 
coupling is broken down because a is large enough to support the velocity perturbation of photons. We checked that 
the Silk damping arises again if we turn off a by hand. We found that the amplitudes of CMB temperature and 
B-mode polarization power spectra depend on the initial amplitude and the aether parameters as Ci oc Avcf 3 c~ 4 . By 
comparing the TT CMB power spectrum induced from the aether field with latest observations, we obtained a new 
constraint on the aether parameters as Avc\ z c~ A < 2.0 x 10~ 17 . 

In addition, we found the dependence of generated magnetic fields on the initial amplitude and the aether parameters 
as \J (B 2 ) cx _4yc 2 3 c~ 8 , which is slightly different from Cgs. We then estimated the amplitude of magnetic fields with 
the aether parameters within the limits of current observations. We found that the maximum amount of magnetic 
fields can be as large as \J (B 2 ) ~ 10~ 22 Gauss at k = 0.1 /iMpc -1 . 

The shape of the magnetic field spectrum is also different from the one with vector perturbations in General 
Relativity 25]. First, the exponential cut-off in the magnetic field spectrum found in [25[ does not arise. The reason 
is exactly the same as in the case for the velocity perturbation of photons mentioned above since the source term, i.e., 
velocities, can survive the Silk damping effect. Secondly, since the vector-mode metric perturbation a decays inside 
the sound horizon, the velocity difference Sv also decays together with a, once the tight coupling between photons 
and baryons breaks down. Thus the spectrum of magnetic fields has a characteristic peak near the Silk damping scale 
at the recombination epoch. Moreover, around the recombination epoch, magnetic fields at larger scales than the 
peak scale experience a little enhancement due to the breakdown of the tight coupling, although the overall spectrum 
shape does not change very much. 

It would be interesting to consider the generation of magnetic fields in other modified gravity theories or other 
inflation models which can amplify the initial power spectrum. For instance, the Einstein- Aether theory can be 
extended to have a more general kinetic term, so called F{K) gravity (34J- It may also be interesting to constrain the 
aether parameters from other cosmological observations, for instance, the correlations or weak lensing effects of large 
scale structure. We leave these subjects for future works. 
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Appendix A: Parameter constraints 



In this appendix, we summarize the current observational and theoretical constraints on the aether parameters. 



1. effective gravitational constant 

To have positive effective gravitational constants, G cos and Gn, 



a < 2 (Al) 
c M > -2 . (A2) 



2. BBN constraint 



The primordial helium abundance created by the Big-Bang Nucleosynthesis (BBN) is affected through the expansion 
rate of the universe via the Friedmann equation. Since the Friedmann equation (Eq. (|16jl) contains the aether 
parameters, the measurements of the primordial helium abundance can constrain the aether parameters as [43| 

ci4 + a < 0.2 . (A3) 



This inequality can be easily satisfied if we set C14 = —a. We assume this relation in our numerical calculations. 



3. PPN limits 



Parametrized post- Newtonian (PPN) parameters in the Einstein- Aether gravity have been derived and constrained 
in refs. |52H54| . PPN parameters in the Einstein- Aether gravity are reduced to two parameters a\ and ai which are 
related to the aether parameters ci, C2, C3 and C4. The constraints from solar-system tests are written in the form as 



~^ + c ^) < 1.7 x 10 -< 
2ci - c\ + C3 " 



ai 

(A4) 



^ ^ a, _ (2 Cl 3-c 1 4)(a + c 14 ) < L2 x 1Q _ 7 

2 Ci23(2-Ci 4 ) 



4. Stability of Scalar- Vector- Tensor perturbations 



The linear perturbations in the Einstein- Aether gravity have been analyzed [42|. Scalar-, vector- and tensor- 
perturbations must be stable in the context of quantum and classical treatments. The classical stability imposes the 
condition that the square of the sound speed must be positive. The quantum stability imposes the condition that 
the ghost should not appear. In other words, the coefficient of a kinetic term of each mode must be positive. These 
conditions imply that 



1 + C13 > , 

CM < 0, Cl < /f 3 , 

2(1 + ci 3 ) 

- 2 < C U < 0, C123 < . 



(A5) 



5. Cherenkov radiation and Superluminal motion 



If the sound speed of metric perturbations is smaller than the speed of light, the transverse-traceless graviton has 
sub-luminal dispersion 55]. Accordingly, relativistic particles will lose its energy by emitting gravitons through a 
similar process of Cherenkov radiation. Therefore we can constrain the aether parameters from the fact that high 
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energy cosmic rays have been observed on the Earth. However, we can avoid these constraints if we assume that all 



modes propagate super-luminally. The conditions of superluminal propagation are given in ref. 42j as 

(2 + ci 4 )ci2 3 < (2 - a)(l + ci 3 )ci4 , 

2c 4 > — ^- , (A6) 

1 + Ci 3 
C13 < . 

The connection of a superluminal propagation and a violation of causality is not trivial and has been still a matter of 
debate. In this paper, we assume that all modes propagate superluminally for simplicity. 

6. Anisotropic stress 

Anisotropic stress of long- wavelength adiabatic modes in scalar perturbations should not be too large [42[ . From 
this we have a constraint as 

|ci 3 |<l. (A7) 

7. Isocurvature mode 

To prohibit isocurvature modes from being dominant at superhorizon scales, we find a constraint (42|: 

— > -1 , (A8) 

Ci4 



which is satisfied if we assume C14 = —a. 



8. Radiation damping 



It is well known that the rate of orbital decay of the Hulse- Taylor binary B1913+16 matches the one induced by 
the emission of gravitational waves in General Relativity (5(| [57j ■ In the Einstein- Aether gravity, the deviation of the 
rate from General Relativity is controlled by the parameter A 



A 



C14 
2 



J_ 2ci4Ci 3 I C14 f 2a 2 - ai \ 2 J_ 

c t (2ci + c\ - elf c v 6(2 + ci 4 )V 2(2ci 3 -ci 4 )/ c s 



(A9) 



where ai and CH2 are the PPN parameters, and c Sjt , it are the sound speed of scalar-, vector- and tensor- 
perturbations 58] . To match the observed rate, the parameter A is constrained as 

\A - l\ < O(10- 3 ) . (A10) 

Appendix B: Cosmological perturbations in Einstein- Aether gravity 

Recently, perturbed Einstein- Aether gravity is discussed and formulated by Refs. [13, S, H3, 5H Ell . Here we 
summarize the cosmological perturbation theory in the Einstein- Aether gravity. 

1. Perturbed equations 

Here, we summarize the perturbed action and equations up to first order in real and Fourier spaces. In the 
synchronous gauge, the metric is given by 

ds 2 = a 2 \-dr] 2 + (<% + hij) cfeW] . (Bl) 



16 



The aether field is written as 



A" 



1 V 1 
a ' a 



(B2) 



Afj, = (-a,aVi) 



First, the equation of motion for the aether field, i.e., Eq. (|6]) leads to following equations. The /Lt = component 
reads 



- a 2 6X = (c 2 + c 3 )V k k - (2 Cl + c 2 + c 3 )UV k k + c 2 ^h k k - (2 Cl + c 2 + 2c 3 )~?tt fe fc , 



(B3) 



where (5A is the variation of the Lagrange-multiplier. We use Eq. (|B3|) to remove <5A from perturbed equations. The 
jtt = i components are 



ci4 + 2ftV$) - (a - ci 4 )HK t + (a + c lA )U 2 V t 



ci^ ' k k + (c 2 + c 3 )V^ fe + c 2 ±h k k>i + ( Cl + c 3 )~A\ fe 



= 



Secondly, we rewrite the Einstein equation (J8j) as 

1 



I? 



-Rg^v = 8nG 



\8irG y 



Then the perturbed Einstein equations lead to (59[: 
h\ + Hh\ = 8iTGa 2 (ST° Q - 5T\) , 

h i j + mh\ + Hh k h S i J - (S*a,-fe* fc + <9 fc <9 fe /i l 3 - <9 fc <9^\. - flj^fc* ) = 167rGa 2 [ST) - i<5\ 

ZHh't + 0- <),!,■ , - d k d k h\ = -\6irGa 2 5T° , 
d 3 h l 3 - tftfj = 16TrGa 2 8T l . 

The perturbed energy momentum tensor for the aether field can be written in the form: 



a 2 5T (A)0 Q = 



c u V\ - (a - c 14 )HV k k - a\uh k k 



+(a + c 14 )H 2 ^ + i( Cl - os) (V fe <* fe - V-' fc fe ) 



C2 (V% + -h\j + 2c 2 [uv k k + -m\ 



We obey the convention of the Fourier transformation as 

d 3 k 



= / ^/(M)e** 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



Furthermore we decompose these Fourier modes into scalar-, vector- and tensor-modes [59(. The definitions of the 
scalar, vector and tensor decompositions are 



A=±l 



(A) 



(B9) 



A=±l 



A=±2 
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where o| A) and 0\f are the projection operators for scalar (A = 0), vector (A = ±1) and tensor (A = ±2) modes. In 
the same way, the energy momentum tensor should be decomposed as 



5T%(k,ri) = 8T^O (0)l + <*T„ (A) (A)< , 

A=±l 

A=±l A=±2 

The perturbed energy momentum tensor for ordinary matter is expressed as 



(BIO) 



6T° 






ST iso 


= -3pn T , 






= - (p + p)v {s) = 




6T t {0) 


= pn^ , 






= -(p+p)« W = 






= pir (v) , 




5T (±2) 


= pttW . 





(Bll) 



Then the Scalar-vector-tensor decomposition gives the following equations. 

Scalar-mode 

Equations of motion for the aether field are 

a 2 S\ = (c 2 + c 3 )kV^ - (2ci + c 2 + c 3 )HkV^ + ^c 2 h iso - ^(2ci + c 2 + 2c 3 )Hh iso , 

akhi so = . 

The scalar components of the energy momentum tensor for the aether field are 
a 2 § Ti so = -a WW +2HkVW + ^h iS o + Hh iso , 

a 2 6T% = - c lA kV {0) + (a- c u )HkV {0) + ^a%h iso , 
a 2 5T^ = - Cl4 (yW + ZHVM) +(a- c lA )HV^ - (a + c u )U 2 V^ , 
a 2 ST t {0) = - l Cl3 + 2Wft(°>) - c 13 fc + 2W(°>) . 

The Einstein equations are 

(l - |) (h iso + Uh^j = (a + c 14 ) (kV^ + UkV^ + 87rGa 2 (p8 + 3pU T ) , 

(l - |) ^Aiso + y (h iso - fc(°>) = -c 14 fcV (0) + (a - c 14 )mV (0) - 8^Ga 2 p<5 , 
fc (% S o - /* (0) ) - -3 [ci4 (^ (0) + 2WV<°>) - (a - Ci 4 )W (0) + (a + c 14 )H 2 y(°)] - 24^GaY s > , 
(1 + c 13 ) (KM + 2Hh^) + y (h iso - ft(°>) = -2c 13 fc (l/W + 2W<°>) + 167rGaV( s ) . 



(B12) 



(B13) 



(B14) 



Vector-mode 

Equation of motion for the aether field is 



ci4 (V^ 1 ) + 2W< ±1 )) - (a - c 14 )WF (±1) + (a + c lA )H 2 V^ + dfcV* 1 ) + \c lz kh^ = . (B15) 
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The vector components of the energy momentum tensor for the aether field are 

a 2 ST (±i) = _ Cl4 (>(±D + 2HV {± V) + (a - c 14 )H7 (±1) - (a + c u )H 2 V^ - \(d ~ c 3 )fc 2 ^ (±1) 
= l Cl3fe 2y(±l) + I Cl3 fcfc(±D , 



(B16) 



and 



a 2 <5T t (±1) 



1 



jCl3 [fci> (±1) + 2^y( ±1 » + ^ + mti^ 

where we used Eq. (|B15|) to derive the second line from the first line in Eq. (IB16[) . The Einstein equations are 

& + 2Ha 



(B17) 



Cl3 <V + ZHv) 



1 + C13 V (1 + Ci 3 )fc 

k 2 a = — - — \l&TrGa 2 q - c 13 k 2 V] , 
l + cia L J 

where we defined new variable a = h^ 1 ' /k. 

Tensor-mode 

The tensor component of the the energy momentum tensor for the aether field is 



The Einstein equation is 



1 + C13 1 + C13 



(B18) 



(B19) 



(B20) 



2. Perturbed action and initial power spectrum 

In this subsection, we calculate perturbed action up to second order to obtain the initial power spectrum of the 
aether field at the end of inflation. For simplicity, the action is decomposed into each component as 



S = Sg + Si + Sa , 



(B21) 



where 



Sn = 



d x^J—gR , 



Sr 



J d i x y /^gCi , 



(B22) 



Sa = 



16ttG 



d^x\J —gC,A 



Here Ci is the Lagrangian density for the inflaton field. Then the perturbed action up to second order is given by 



:(2) 



1 



16ttG 



rf 4 :; 



uk It 



h kl h M - h k k h e e + [m 2 + 57i) h k k h\ + (-m 2 - m\ h M h kl 



(B23) 



S) 



(2) 



. d 4 x- 
16ttG / 4 



j [(2H 2 + m) {h kt h M - h k k h\) -a(H 2 + ZH) (h M h M - h k k h\) 



(B24) 
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?(2) _ 



16ttG 



c? 4 xa 2 



-c u V*Vi + (a + c 14 )U 2 V l V, - (a - ci 4 )WV*V5 



1 



a ( U z + 2H ) ft M /i. 



n [V t ' k + ~h t k 



I;, 



I'M 



■a[U z - HJ h\h l t 

h\)+ C2 (v. + h' 



V 



S I h k 



(B25) 



+c 3 V, + V\ + —h 



1:. 



Here, we employed the single field slow-roll inflation model and replaced the inflaton field with T-L and metric pertur- 
bations by using the Eqinstein equations and the equation of motion for the inflaton. 

Hereafter we focus on the vector-mode only. By performing Fourier transformation and scalar- vector-tensor decom- 
position, we have 

d 3 k 



—eH 2 \v\ 2 -c 2 v k 2 \v\ 2 

C14 



(B26) 



where v = y/—ci4aV/\/8irG and Ti 2 —H = eH 2 . Through variation with respect to v, we have the equation of motion 
for the vector perturbation as 



v k + c%k l v k + — eH z v k = 

C14 



(B27) 



Supposing that the slow-roll parameter e is constant during inflation, the conformal time 77 and the scale factor a are 
expressed as 



n 



1 1 a 
H 1 — e' ai 



1/(6-1) 



(B28) 



where 771 and ai are the conformal time and scale factor at the end of inflation. Then we can solve Eq. (|B27|) easily 
to obtain 



where V[ n { 



v k (v) = JjC-^/V^+^/^w (-fav?) , 



and is the Hankel function of the first kind. Finally the power spectrum of the 



(B29) 



4 C14 (1— e) 

aether field at the end of inflation is given by 



2ir 2 



(V(r,,-k)V*( V X)) = (27rf- T Vv(k)S(k-k') 



k 3 



where 



Tv(k,m) 



{l-ef 



^r'M (f) (-W 



(B30) 



(B31) 



-C14 (8ttG)- 

with kg = 0.002Mpc _1 and n v — 3 — 2v m i. From the above equation, we can read off the amplitude of the initial 
power spectrum as 



A v = 



{l-e) 2 H{ 



^ M (|)- 2 ^(-^r 



-c l4 (8ttG) 



(B32) 



If we substitute the above amplitude into Eq. (|49|) with the condition that C14 = —a, the inequality Eq. (|49|) can be 
rewritten as 

^2 



< 10 



4:i 



(B33) 



where we have assumed e ~ 0.016, which corresponds to the scalar-tensor ratio r ~ 0.01, H\ ~ 5 x 10 13 GeV, and 
m ~ -2.0 x 10 15 GeV. 
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